Nuclear structure corrections in muonic deuterium 
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The muonic hydrogen experiment measuring the 2P — 25 transition energy [R. Pohl et al., Nature 
466, 213 (2010)] is significantly discrepant with theoretical predictions based on quantum electro- 
dynamics. A possible approach to resolve this conundrum is to compare experimental values with 
theoretical predictions in another system, muonic deuterium /iD. The only correction which might 
be questioned in is that due to the deuteron polarizability. We investigate this effect in detail 
and observe cancellation with the elastic contribution. The total value obtained for the deuteron 
structure correction in the 2P — 2S transition is 1.680(16) meV. 
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The determination of electromagnetic properties of nu- 
clei from precise atomic spectroscopy has become possi- 
ble due to significant progress in atomic structure the- 
ory, which in turn is based on quantum electrodynamics 
(QED). The proton charge radius as obtained from the 
Lamb shift in hydrogen [l| is more accurate than any de- 
termination using electron scattering, including the most 
recent ones 0, [j], and in agreement vi^ith them. The 
determination of the deuteron charge radius from the 
measurement of the hydrogen-deuterium isotope shift in 
the 2S-1S transition frequency d, Q, apart from being 
the most accurate, has stimulated reanalysis of the elec- 
tron scattering data. At present the atomic isotope shift 
determinations of charge radii for helium, lithium and 
beryllium atoms are not only the most accurate ones, 
but also the only ones available for short lived isotopes 
i. 

It was a great surprise that the proton charge radius 
Tp determined from the muonic hydrogen Lamb shift [3] 
gave a result in conflict with the value determined from 
electronic hydrogen. The 5 a discrepancy in rp is the 
first indication that our knowledge of interactions in these 
simple atomic systems is not complete. While we will not 
pursue the possible explanations of the proton charge ra- 
dius discrepancy, we point attention to the critical test 
which can be performed with muonic deuterium. The 
electronic H-D isotope shift gives a very accurate differ- 
ence for the deuteron and proton charge radii. If the 
results for the difference between /iD and are consis- 
tent with the electronic H-D, this would mean that there 
is an extra muon-proton interaction, which cancels out in 
the /iII-/iD difference. In order to draw these conclusions, 
all other effects contributing to the isotope shift have to 
be analyzed. The only correction which goes beyond the 
standard QED treatment is that due to the nuclear polar- 
izability, and is in general due to the nuclear structure. 
It is the purpose of this work to study these effects in 
muonic deuterium. 

Before this, however, we note that there is no unique 
definition of the charge radius for all nuclei. It depends on 
the nuclear spin, and in particular, there is no established 
definition for spin 1 nuclei, such as the deuteron. The 
mean square charge radius (-R^) for an arbitrary spin / 



particle is defined through the effective interaction with 
the electromagnetic field 
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where fi and Q are the magnetic dipole and the electric 
quadrupole moments. For a scalar particle (5o = 0, and 
for a half-spin particle S1/2 = 1/8. For a vector and 
higher spin particle we proceed as follows. The most 
general Lagrangian for particle with spin / = 1, which 
includes terms linear in the electromagnetic field strength 
FisM 
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where u^^, = 
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pt^u ^ dt^A^-d^At". The effective nonrelativjstic Hamil- 
tonian obtained from Eq. ^ gives the V • E term with 
the coefficient eQ/6 [9|] . What part of this coefficient 
should be included in the charge radius, and what in the 
kinematical term Si ? The most natural assumption is 
that {R'^) = for a point 1 = 1 particle. One possible 
choice for a point vector particle is 5 = 1, Q = 0, and an- 
other one g = 2, Q = —l/rn?. The second choice has the 
advantage that it leads to the renormalized QED theory 
of a charged vector particle, while the first choice leads 
to the simplest form of the Lagrangian. In this work, fol- 
lowing [lOj, we adopt the first choice, and consequently 
assume for a vector particle 5/ = 0. This assumption 
affects the relativistic recoil correction, while the finite 
nuclear size correction is 
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where (/>^(0) = {rrir of / [n n?) 5iq, and is the /iD re- 
duced mass. 
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Having defined the leading finite size effect, we pro- 
ceed to tlie evaluation of nuclear structure corrections. 
A nucleus is not a rigid particle, it can be excited by an 
orbiting electron or muon, which results in the shift of 
atomic energies. In the following, we derive general for- 
mulas for the nuclear polarizability shift with any Hamil- 
tonian for deuterium, using the perturbation expansion 
in the muon mass m over the deuteron mass niD, includ- 
ing the so called Coulomb correction. Here, the deuteron 
binding energy counts as vn? /iriD. With the assumed 
accuracy relativistic effects for the muon, as well as for 
the deuteron, can be treated perturbatively and we start 
derivation from the leading electric dipole excitations. 
The nonrelativistic formula for the electric dipole nuclear 
(scalar) polarizability correction is 



where «£; is the electric dipole polarizability of a nucleus 

and E is the mean excitation energy. Both of them 
have already been accurately calculated for the deuteron, 
namely = 0.6330(13) fm^ and E = 4.94 MeV fll|. 
Our numerical results obtained in this work are in an 
agreement with them. We shall mention that the for- 
mula for the leading Coulomb correction in 15* and 25* 
states was first obtained by Friar in [T2l |. 

The next to leading Coulomb correction, which has 
not been considered so far, written explicitly for 2P — 2S 
transition, is 
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where Hq is the nonrelativistic Coulomb Hamiltonian for 
the muon with reduced mass to,.. Denoting the nuclear 
excitation energy by E, the polarizability correction is 
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The nuclear excitation energy E is much larger than 
a typical atomic (muonic deuterium) excitation energy, 
thus one may perform the large E expansion of the elec- 
tronic matrix element. The appropriate formula for this 
expansion in atomic units, £ ~ E/(rnr a^) is 
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From this expansion, the leading electric dipole polariz- 
ability correction is |13 | 
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while the leading Coulomb correction, written explicitly 
for 2P — 25* transition, is 
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The dipole operator d is the position R of the proton with 
respect to the deuteron mass center. It is thus assumed 
that there are no corrections to the electric dipole opera- 
tor, and for example mass center of the proton coincides 
with the charge center of the proton within the nucleus. 
The uncertainty introduced by this approximation is un- 
known. This is due to the fact that the underlying QCD 
theory is nonperturbative and the resulting exact theory 
of nuclear forces is not yet known. 

In the evaluation of further corrections we take the in- 
finite nuclear mass limit, thus neglecting nuclear recoil. 
The relativistic corrections to the electric dipole polariz- 
ability effects can be obtained from the two-photon ex- 
change amplitude [l3j 
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The nuclear excitation energy E is much smaller than the 
muon mass to and the nonrelativistic contribution comes 
from the region lo ^ E and k « \j2mE. Thus to obtain 
it one neglects uP' and w^. The leading nonrelativistic 
term agrees with that in Eq. (|8]), with the reduced mass 
TOr replaced by a muon mass to. Relativistic correction 
comes from the next terms in the small E expansion, 
namely 

5bE = ze4(/.2(0)i / dE{<t,n\d\E? f 
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This is the only relativistic correction which is not negh- 
gible at our level of accuracy. 

The corrections due to higher multipole polarizabilities 
and higher order corrections due to the finite deuteron 
size can be treated together, and we show that they can- 
cel each other at the leading order. Let us consider at 
first the related electronic matrix element P for the non- 
relativistic two-photon exchange 
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where is the nonrelativistic Hamiltonian for the muon 
in the nonrecoil limit, and i? is a position of proton with 
respect to the nuclear mass center. Using on mass shell 
approximation with subtractions of the leading Coulomb 
interaction, the finite size, and the electric dipole contri- 
bution, it becomes 
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where we performed expansion in the small parameter 
^/2mE\R — R'\. The corresponding correction to atomic 
energy is 

hE = - ^dEjd'Rd'R' (f^UR) <t>E{R) <t>D{R') c^UR) P 

(16) 

Consider now the first i?-independent term. When 0_e = 
(t>D, it will be the elastic part know as a Zemach correc- 
tion [l^, but the inclusion of all excited states leads to 
5{R— R') and the \R — R'\^ term vanishes completely. As 
a result, there is no Zemach [elastic 0{a^)] correction for 
muonic deuterium. From Eq. (jTS]) only the second term 
remains, which gives 
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These corrections are due to the electric dipole, 
the quadrupole and the monopole nuclear excitations, 
namely 
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As we assumed at the beginning, all corrections of or- 
der TV? jmY) are neglected. However, due to the large 
magnetic moment anomaly of the proton and the neu- 
tron, we make an exception and consider the magnetic 
dipole polarizability correction. It comes from 
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where gp = 5.586, and g„ — —3.826. We use the anal- 
ogous two-photon exchange formula as for the electric 
dipole transitions 
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and perform the nonrelativistic approximation 
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Corrections due to the intrinsic proton polarizability and 
the proton Zemach moment are to a good approximation 
the same as in muonic hydrogen. Therefore, we use a 
recent result of very thorough calculations IS] AE{2S) = 
—36.9(2.4) /ieV and scale it by a factor 
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The final expression for the nuclear polarizability com- 
bined with the elastic contribution is 

AE = 5oE + 5cE + SpE + SqE + 6mE + 5pE. (23) 

Numerical calculation of deuteron matrix elements 
are performed using the discrete variable representation 
(DVR) [III method. In the DVR method, the Hamilto- 
nian is represented as a symmetric matrix, which can be 
diagonalized and all formulas represented as a finite sum 
over the spectrum. Numerical results using the mod- 
ern nucleon-nucleon AV18 potential from Argonne Na- 
tional Laboratory [17,] , are presented in Table I. We have 
checked numerics by the calculation of the electric dipole 
polarizability, and out result ad = 0.634 fm'^ is close 
to the recommended value from [ll|. The difference of 
about 0.001 fm'^ comes from the fact that we take into 
account a small neutron-proton mass difference in the 
electric dipole operator. 

Surprisingly, the total value AE is close to the nonrel- 
ativistic electric dipole polarizability contribution SqE + 
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TABLE I: Nuclear structure corrections in muonic deuterium 
for 2P — 2S transition. 



correction 


value in mcV 


SoE 


1.910 


SciE 


-0.255 


Sc2E 


-0.006 


5rE 


-0.035 


s: IT' 

oqoE 


—0.045 


SqiE 


0.151 


5q2E 


-0.066 


SmE 


-0.016 


5pE 


0.043(3) 


AE 


1.680(16) 



dciE + SciE — 1.649 meV. This means that relativistic 
and higher multipole corrections, although individually 
not small, tend to cancel between themselves. For the 
final uncertainty, we assume 50% of these higher order 
corrections. We can not at this moment give a more reli- 
able estimate of uncertainty, but note that it is about 20 
times smaller than the discrepancy in muonic hydrogen. 

Our result for the nuclear structure correction AE is 
not in good agreement with former calculations. Lei- 
demann and Rosenfelder in [l^ obtained for the po- 
larizability correction of 2S state the result —1.500(25) 
meV. This should be combined with the elastic contri- 
bution obtained by Martynenko [l9] of —0.37 meV and 
the Coulomb correction of 0.26 meV, which totals to 
— 1.61(3) meV. This is 2tT away from our result, shown 
in Table I. The difference may come from three sources. 
The first one is lack of clear separation between the elas- 
tic contribution from [l^ and the inelastic one from [l^ , 
thus some terms might be counted twice. The second one 



is neglect in [iSj of the intrinsic proton polarizability cor- 
rection of 0.013 meV [l^. The third source is the extra 
coefficient i?^ = 0.9778 used in [l^ for the polarizabil- 
ity correction, which reflects the fact that probability of 
finding muon within the deuteron is not exactly (f)'^ (0) but 
Rf^ (/'^(O). To verify this coefficient one has to investigate 
three photon exchange correction, details are beyond the 
scope of this work but we claim lack of such coefficient. 

In order to shed light on the proton charge radius 
discrepancy, we consider the difference Ed{2P — 2S) — 
Eh{2P — 2S) m^jj/mfjj, where the proton size and the 
proton polarizability cancel out. This difference is sensi- 
tive to the deuteron structure radius, which is known 
from a very accurate H-D(2S' — IS) isotope shift. If 
agreement between experiment and theoretical predic- 
tions based on QED calculations [13, HH including nu- 
clear polarizability correction calculated in this work is 
observed, this may mean that the muonic hydrogen dis- 
crepancy is caused by a local (^ fm) muon-proton in- 
teraction or by a 5cr shift in the Rydberg constants [7|. 
There could be as well a different source of the discrep- 
ancy, a long-range type interaction (^ 200 fm), and this 
will not cancel out in the difference (??), and a small 
discrepancy in /xD will persist in this case. 

In summary, we have demonstrated in this work that 
accurate predictions for muonic deuterium are feasible, 
the non-QED corrections have been here accurately cal- 
culated, and comparison with experimental transitions in 
fiD will give hints on a possible source of discrepancies 
in the proton charge radius. 
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